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rule. 



Radiative decays of the radially excited charmonium res- 
onance, if)(2S), into 7T7T, KK and r/rj final states have been 
measured in a sample of 3.96 x 10 6 tp(2S) events collected 
by the BES collaboration. The branching ratios B(ip(2S) — > 
7/ 2 (1270)) = (2.27 ± 0.26 ± 0.39) x 10" 4 and B(ip(2S) -> 
7/o(1710)) x B(/o(1710) -> K+K') = (5.59 ±1.12 ±0.93) x 
10" 5 are obtained. When compared to the corresponding ra- 
diative J/tp decays, the observed ip{2S) radiative decay rates 
into 7/2(1270) and 7/0 (1710) are consistent with the "15%" 



I. MOTIVATION 

In perturbative QCD, the dominant process of J/tp and 
ijj(2S) hadronic decay is cc annihilation into three gluons. 
Since the decay width is proportional to the amplitude of 
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the cc wave function at the origin, ^(O)! 2 , the branching 
fractions of J/ip and ip(2S) decays into light quark states 
are related as : 

B(i>(2S) -+h) B(iP(2S) -> ggg) 



B(J/ip -> ft) 
B(ip(2S) -> e+e- 
B(J/tP ^ e+e^ 



= (14.6 ±2.2)% 



This relation is called i/ie i<5% ra/e. The prediction 
was originally made for the total decay width into three 
gluons. Since the partial widths of individual channels 
involving the initial annihilation of cc quarks arc also 
functions of the |^(0)| 2 , we expect this rule to be gener- 
ally valid. 

Results from the Mark II experiment pi show that 
while many of the ip{2S) hadronic decay channels obey 
this rule, it is severely violated in vector plus pseudo- 
scalar (VP) final states such as pn and K*K — the 
so called pm puzzle. The BES experiment also reported 
heavy suppression in the vector plus tensor (VT) final 
states such as K*K%, pa-i, ajf 2 and </>/ 2 [f|. 

In perturbative QCD, the radiative J/tp and ip(2S) de- 
cays should be similar to hadronic decays except instead 
of decaying into three gluons, the radiative mode decays 
via two gluons and one photon. Thus one power of the 
coefficient as is replaced by cxqed in the cross section 
formula. It is expected that the "15%" rule should also 
work for radiative decay modes M. Hence the ratio of 
B(ip(2S) -> jX) to B(J/ip -> jX) for different final 
states X should be roughly 15%. This paper explores 
the ij)(2S) radiative decays into pairs of pseudoscalars, 
7T7T, KK and rjrj, and reports the first branching frac- 
tion measurements of ip(2S) — > 7/2(1270) and 7/o(1710). 
The branching fractions of Xco an d Xci decays into 7T7r 
and 7777 are also reported. 

The /o(1710) has been observed with a large branching 
ratio in radiative decays of J/tp into KK, but not in the 
reaction K~p — > KKA by the LASS experiment ||. The 
later result excludes /o(1710) as a conventional ss state 
and makes it a leading glueball candidate. Thus whether 
/o(1710) can be seen in the radiative ip(2S) decays is 
quite interesting. 



II. BES DETECTOR 

This study uses a subset of the 3.96 million e + e~ — > 
ip{2S) event jq] logged by the BES detector operating 
at the BEPC storage ring. A detailed description of 
the BES detector can be found elsewhere Jl|]. It fea- 
tures a 40-layer main drift chamber (MDC) in a 0.4 T 
solenoidal magnetic field providing a momentum reso- 
lution of cTp/p = 1.7% yl +p 2 (GeV/c) and a dE/dx 
resolution of 9% for hadron tracks. Outside of the 
MDC cylinder is an array of 48 scintillation counters 



of the time-of-flight (TOF) system with a resolution of 
450 ps for hadrons. A 24-layer lead-gas barrel elec- 
tromagnetic shower calorimeter (BSC), outside of the 
TOF system, provides an energy resolution of cte/E = 
0.22/ y/E (GeV), and spatial resolutions of 0$ — 4.5 mrad 
and ug = 12 mrad. The BSC is surrounded by a magnetic 
coil and steel plates (for magnetic flux return). There is 
a 3-layer p counter interleaved with the steel plates to 
identify p tracks. 

The photons used in this analysis are detected as show- 
ers in the BSC. Showers within 8° of each other are 
regarded as split showers of a single photon candidate 
and their energies are recombined. For a 7r° or 77 de- 
caying into two photons, the helicity angle of the de- 
cay should be flat and hence | cos dhelitity | < 0.99 is re- 
quired to remove some asymmetric background decays. 
For charged tracks, only those that fall in the fiducial 
region | cos f?| < 0.8 are used. Charged tracks identi- 
fied as electrons or positrons by the BSC or identified as 
p + or pr by the p counter arc rejected. TOF informa- 
tion, dE/dx information and kinematic fitting are used 
to identify charged particles. Kinematic fits are applied 
to improve momentum measurements and mass resolu- 
tions, as well as to resolve combinatorial ambiguities if 
more than one combination is possible in the same event. 



III. ip(2S) -»■ 77T7T ANALYSIS 

Here the analysis of decays into both charged and neu- 
tral pion pairs, tp(2S) — > 77r + 7r~ and ip(2S) — > 77r°7r°, 
is described. In the charged channel, events are se- 
lected with two oppositely charged tracks and at least 
one photon. Kinematic fit, TOF, and dE/dx informa- 
tion are combined and the probability for the two pions 
hypothesis should be greater than 1% and should also be 
greater than the probability for the two kaons hypoth- 
esis. The solid-line histogram of Fig. [l] is the ir + ir~ 
invariant mass distribution below 2.5 GeV. A clear p sig- 
nal is observed over a continuous background. These are 
due to the initial state radiation processes e + e~ — * 7/5 
and e + e~ — > "fp + p~. These processes are also pre- 
sented in the e + e~ — > tt scan data taken by BES at 
v/s = 3.55 - 3.6 GeV (below ip(2S) threshold) (To) and 
the mass distribution from the tt scan can be used to rep- 
resent the background in the ip(2S) sample. The n + ir~ 
mass distribution taken from the r scan data is shown 
in Fig. [j] (with the dashed line). The backgrounds are 
removed by assigning each ip(2S) event a weight of 1 and 
each r event a weight of —w in the likelihood function 
for the mass fit. Here w is the ratio of the integrated 
luminosities of these two data sets. Fig. ||a shows the 
result of subtracting the tt scan histogram normalized 
to the ip{2S) from the ip(2S) histogram. 
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M(7T + 7T ) 

FIG. 1. M^ +7r - from ip(2S) data in the solid line his- 
togram and luminosity normalized r data in the dashed line 
histogram. 

In the neutral mode, events are required to have at 
least 5 neutral tracks and no charged tracks. A six- 
constraint fit is made to all possible 77r°7r° combinations 
with two 7T° resonances. The combination with the small- 
est fit x 2 is selected. A four-constraint fit is also applied 
on that combination and |M 77 — M^o\ < 70MeV is re- 
quired. The resulting mass distribution is shown in Fig- 
ure |^b. 



Branching fractions of B(^(2S) -> 7/ 2 (1270)) = 
(2.21 ± 0.26 ± 0.39) x 10~ 4 from the ip(2S) -> ~/ir + ir- 
channcl and B(^(2S) -> 7/ 2 (1270)) = (2.95 ± 1.10 ± 
1.12) x 10" 4 from the ip(2S) -> 77r°7r° channel are ob- 
tained using the total number oftp(2S) events and detec- 
tion efficiencies, which are determined from Monte Carlo 
simulations, of 43.8% and 9.6%, respectively. The com- 
bined result from these two channels is 

B(4>(2S) -> 7/2(1270)) = (2.27 ± 0.26 ± 0.39) x 10~ 4 . 

Here and below, the first error is statistical and the sec- 
ond is systematic. The latter includes uncertainties from 
varying the cuts, the shape of the background, and the 
uncertainty from the total number of ip(2S) 's. This result 
is consistent with the "15%" rule when compared with 
the corresponding branching fraction from J/ip decay on 
PDG. See Table |. 

A /o(1710) signal is observed in the 7r + 7r~ invariant 
mass distribution, as shown in Fig. ||a. The number of 
/o(1710) events above background is 39.2 ± 10.1. The 
Monte Carlo determined efficiency for this channel is 
45.4%, and the resulting branching fractions are 

B(il>{2S) -> 7/o(1710)) x fl(/ (1710) -> tttt) 

= (3.38 ± 0.87 ± 1.41) x 10" 5 




l(nV) 



FIG. 2. (a): M v +„- fit result. The data points are ob- 
tained from the difference of the two histograms in Fig. [l] (b): 
M n o n o. 

In both the ir + ir~ and n ir invariant mass distribu- 
tions, clear /2(1270) signals are observed. Both distri- 
butions are fitted with a D-wave Breit-Wigner function 
with the resonant parameters fixed at the PDG || val- 
ues for the /2(1270). An S-wave Breit-Wigner with mass 
and width fixed at the PDG values for the /o(1710) is 
included in the mass fitting in the 77r + 7r~ channel in or- 
der to describe the line shape in that region. A phase 
space background is included in the 77r°7r° channel. The 
fit yields 209.8±24.9 events and 29.9± 11.1 events above 
background as shown in Fig. ||a and Fig. ||b, respectively. 



(90% C.L. 0) 



< 5.50 x 10~ 5 



The region with a 7r + 7r~ invariant mass greater than 3 
GeV has been presented elsewhere |Q . The region with a 
7r°7r° invariant mass greater than 3 GeV (see Fig. ||a) has 
signal peaks due to the XcO and % c2 charmonium states. 
This mass distribution is fitted with two Breit-Wigner 
resonancess plus a polynomial background function, and 
96.9 ±11.1 and 20.8 ±5.8 events are obtained for XcO an d 
Xc2, respectively. The detection efficiencies are 10.5% 
and 8.2%, and the resulting branching fractions are 




FIG. 3. Invariant mass of (a) 7r°7r° and (b) r\r). 
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B(xco -> 7r u 7r u ) = (2.65 ± 0.30 ± 0.58) x 10" 



and 



B( X c2 - ^TT ) 



.7 ±2.4 ±5.0) x 10" 



The detection efficiencies, branching fraction accep- 
tances, and final results for the decays described in this 
section are summarized in Tables [Tl| to |v[ 



IV. ip(2S) -► -yKK ANALYSIS 

Here the analysis of decays into charged and neutral 
kaon pairs, tp(2S) -> ~/K+K~ and ip(2S) -> -yK s K s -> 
747r ± , is presented. For the i[>(2S) — > jK + K~ channel, 
cuts similar to those for the 77r + 7r~ analysis are used, 
but with the requirement that the probability from par- 
ticle identification for each kaon hypothesis should be 
greater than 0.01 and should be greater than the prob- 
ability for the pion hypothesis. QED backgrounds such 
as e + e~ — * 70 — * jK + K~ and e + e~ — > 7/i + /j _ are de- 
termined using the r scan data (See Fig. |). A / (1710) 
signal and a hint of a possible / 2 (1525) signal (Fig. ||) 
are observed. The mass distribution is fitted using S- 
wave and D-wave Breit-Wigner functions with masses 
and widths fixed at the PDG values for /o(1710) and 
/a(1525), respectively. The fit yields 71.9 ±14.4 / (1710) 
events above background. The detection efficiency is 
33.6%, giving a branching fraction 





M(K + K~) 

FIG. 4. M K + K - from ip(2S) data in the solid line his- 
togram and luminosity normalized r data in the dashed line 
histogram. 

B(iP(2S) -» 7 /o(1710)) x B(/o(1710) -> K+K~) 

= (5.59 ± 1.12 ±0.93) x 10" 5 

The systematic error includes the uncertainties from 
varying the cuts and the total number of ?/>(2S)'s. This 
result is again consistent with the "15%" rule. See Table 

1 



For the ip(2S) -> jK s K s -> j-k+tt-tt+tt- channel, 
events with two positive charged tracks, two negative 
charged tracks and at least one photon are selected. Both 
Kg vertices are reconstructed and \M„+„- — M K o\ must 
be smaller than 20 MeV. At least one of the combina- 
tions must yield a four-constraint fit with \ 2 probability 
greater than 0.01. If more than one combination survives, 
the combination with the smallest value of 

is chosen. An S-wave Breit Wigner plus a polynomial 
background are used to fit the invariant mass distribution 
shown in Fig. ^p. The fit finds 6.8 ±3.1 events, or 
an upper limit of 10.8 events at 90% confidence level. 
The detection efficiency for this channel is 18.0%. The 
branching fraction is 

B(V>(2S) -> 7/o(1710)) x fl(/ (1710) -> KgK g) 

= (2.10 ± 0.96 ± 1.11) x 10~ 5 



< 3.98 x 10~ 5 (90% C.L.) 

The detection efficiencies, branching fraction accep- 
tances, and final results for the decays described in this 
section are summarized in Tables [0] to |yj 
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FIG. 5. Invariant mass of (a) K + K~ and (b) K° S K%. The 
data points in (a) are obtained from the difference of the two 
histograms in Fig. 



V. xp(2S) -»■ 77777 ANALYSIS 

Here the analysis of decays into 77 pairs, i/)(2S) — > 
77/7; — > 57, is presented. The selection criteria for this 
channel are similar to those used in the ip(2S) — > 77r°7r° 
channel except that the requirement |M 77 — M n \ < 70 
M eV is imposed. In the region where M m < 3 GeV, no 
7777 resonant state is observed. In the mass region with 
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M nn > 3 GeV, shown in Fig. ||b, a Xco signal and a possi- 
ble Xc2 signal are observed. Two Breit-Wigner functions 
with a polynomial background are used to fit the mass 
distribution, and 12.7 ± 5.3 events for Xco and an upper 
limit of 5.9 \ci events are found, giving the following 
branching fractions 

B{Xco -> w) = (1-94 ± 0.81 ± 0.59) x 10~ 3 

B(xc2 -» rpn) < 1-22 x 10~ 3 (90% C.L.) 

The detection efficiencies for these two channels are 
11.9% and 10.5%, respectively. 

Flavor SU(3) symmetry predicts that the branching 
fractions of Xco decay into 7r°7r° and r\r\ should be the 
same except for a phase space factor and a barrier factor 
of p( 2s + 1 ) i where p is the momentum of the n° or 77 in x c 's 
rest frame and s is the spin of the Xc- Based on the PDG 
values for the Xco, this predicts B(x c o — * w)/-B(XcO —> 
7r°7r°) = 0.95 which is consistent with our measurement 
of 

B(Xco VV) = J3 ± 32 ± 27 
B( X cO -» 

The detection efficiencies, branching fraction accep- 
tances, and final results for the decays described in this 
section are summarized in Tables [n] to [v]. 

VI. rp(2S) NORMALIZATION 

The total number of ip(2S) events in the BES data 
sample is determined from the observed number of cas- 
cade decays of ip(2S) — > tt + tt~ J ftp, J/ip — ► X. The total 
number of tp{2S) — ► ir + ir~J/ip events observed in the 
full BES data sample is (1.227 ± 0.003 ± 0.017) x 10 6 
§ @. The decay branching fraction of B(ip(2S) 
TT+TT-J/tp) = (31.0 ± 2.8)% taken from PDG 2000 gj 
is used to calculate the total number of produced ip{2S) 
events in the data sample. This paper provides the ratios 
of the branching fractions and B(ip(2S) — > 7r + 7r~ J/ip) 
in order to isolate the systematic error caused by the 
ip{2S) — > J/ip branching fraction. See Tables [il] to 

^ for the results. 

VII. SUMMARY 

This paper studies ip(2S) — > r yn + TT~ , jir°ir°, -fK + K~, 
■yKgKg, -fr/ri final states and reports the first mea- 
surement of the ip(2S) -> 7/2(1270) and ip(2S) -> 
7/o(1710) — ► -fK + K~ and 7X^X5 branching fractions. 
A clear /o(1710) signal in tp(2S) radiative decay into 
K + K~ final states is observed. The results are consistent 
with the "15%" rule. 



In addition, this paper reports the first measurement 
of the branching fractions of Xco and Xc2 decay into ir°ir°, 
XcO decay into rji], and an upper limit of the branching 
fraction of Xc2 decay into rjrj. The results from Xco —> 
7r°7r° and r\r\ are consistent with the prediction by SU(3) 
flavor symmetry. 
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G 



Final state 


B{i>{2S) ->)(xHT 4 ) 


B(J/i/> -)(xl0- 4 ) 


B(tP(2S))/B(J/^) 


7/2 (1270) 


2.27 ±0.26 ±0.39 


13.8 ± 1.4 


(16.4 ±3.1)% 


7/o(1710) -> -yK+K- 


0.559 ±0.112 ±0.093 


4.251^ ft 


(13.2+to)% 



TABLE I. Verification _of the 15% rule. The value of B(J/ip -> 7/0 (1710) -> 7A+IO is obtained from 
B(J/ip -> 7/o(1710) -» 7AA) = 8.6i ;9 divided by a factor of 2 to account for isospm. 



Mode 


Number of events 
from mass fitting 


Detection 
efficiency 


Branching Fractions 
correction factor 


Number of events 
after correction 


ip{2S) -> 7/2(1270) from 77r+7r~ 


209.8 ± 24.9 


43.8% 


0.847 ± 0.024 


565.5 ±67.1 ±85.7 


^(25) -> 7/2(1270) from 77r u 7T u 


29.9 ± 11.1 


9.6% 


0.827 ± 0.028 


376.7 ±139.9 ±138.1 


ip(2S) — > 7,/o(1710) — > 77T7T from 77r + 7r~ 


39.2 ± 10.1 


45.4% 


1.0 ±0.0 


86.3 ±22.2 ±35.2 


^(2S) -> 7/o(1710) -> -yK+K- 


71.9 ± 14.4 


33.6% 


1.0 ±0.0 


214.0 ±42.9 ±30.6 


V(2S) - 7/o(1710) - 


6.8 ±3.1 


18.0% 


0.471 ± 0.004 


80.3 ±36.6 ±41.9 



TABLE II. Numbers of events before and after correction for efficiency and branching fraction acceptance. Branching fraction 
acceptance factors include branching fractions from intermediate decays processes such as ty° — » 77, 77 — > 77, Kg — > ir + ir~ , 
/2(1270) — > 7T7T, etc. They do not include iso-spin factors. 



Mode 


P(xHT 4 ) 


B/B{iP(2S) -> 7r+7r-J/7/>)(xlO- 4 ) 


ip(2S) -» 7/2(1270) from 77T+7T- 


2.21 ±0.26 ±0.39 


7.14 ±0.85 ±1.09 


ip(2S) -> 7/2(1270) from 771- V J 


2.95 ± 1.10 ± 1.12 


9.52 ± 3.53 ± 3.49 


tl)(2S) -> 7/2(1270) from 77T7T 


2.27 ±0.26 ±0.39 


7.31 ±0.83 ±1.04 


ip(2S) — > 7,/o(1710) — > 77T7T from 77r + 7r~ 


0.338 ±0.087 ±0.141 


1.09 ±0.28 ±0.44 


V>(2S) — ► 7/o(1710) -> 7A+a- 


0.559 ±0.112 ±0.093 


1.80 ±0.36 ±0.26 


i'i^S) — > 7,/o(1710) — * 7^5 A'g 


0.210 ±0.096 ±0.111 


0.68 ±0.31 ±0.35 



TABLE III. Branching fractions and ratios of branching fractions (B/B(ip(2S) -> tt + iy~ J/tp)) for V(2S) -> 7A -> 7PP 
modes (P stands for pseudo-scalar). 



Mode 


Number of event 
from mass fitting 


Detection 
efficiency 


Branching Fractions 
correction factor 


Number of events 
after correction 




96.9 ± 11.1 


10.5% 


0.9761 ± 0.0005 


945.4 ± 108.3 ± 163.8 




20.8 ±5.8 


8.2% 


0.9761 ± 0.0005 


259.9 ± 72.5 ± 145.8 


XcO -> W 


12.7 ±5.3 


11.9% 


0.1547 ±0.0035 


689.9 ± 287.9 ± 188.9 


Xc2 -» W 


< 5.9 


10.5% 


0.1547 


< 363.3 



TABLE IV. Numbers of events corrected for efficiency and branching fraction acceptance for Xc decay. 



Mode 


P(xl0 _tl ) 


B x B{i/>(2S) ^7Xc0, 2 )(xl0- 4 ) 


B x B{i>{2S) -» j X cO,2)/B(^(2S) -> 7r+7r-J/7/>)(xlO- 4 ) 


XcO - *V 


2.65 ±0.30 ±0.58 


2.47 ±0.28 ±0.49 


7.96 ±0.91 ±1.38 


Xc2 - ^V J 


0.87 ±0.24 ±0.50 


0.68 ±0.19 ±0.39 


2.19 ±0.61 ±1.23 


Xco -» m 


1.94 ±0.81 ±0.59 


1.80 ±0.75 ±0.52 


5.81 ±2.42 ± 1.59 


Xc2 -> vv 


< 1.22 


< 0.95 


< 3.06 



TABLE V. The Xc decay branching fractions, and ratios of branching fractions for Xco.2 — > it it or 7777 
(B x B(i>(2S) -» 7X c o, 2 )/B(V(25') -» tt+tt" J/^)). 
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